Asymptotically exact trial wave functions for yrast states of rotating Bose gases 

S. Viefers and M. Taillefumier 
Department of Physics, University of Oslo, P.O. Box 1048 Blindern, N-0316 Oslo, Norway 

(Dated: May 17, 2010) 

We revisit the composite fermion (CF) construction of the lowest angular momentum yrast states 
of rotating Bose gases with weak short range interaction. For angular momenta at and below the 
single vortex, L < N, the overlaps between these trial wave functions and the corresponding exact 
solutions increase with increasing system size and appear to approach unity in the thermodynamic 
limit. In the special case L — N, this remarkable behaviour was previously observed numerically. 
Here we present methods to address this point analytically, and find strongly suggestive evidence 
in favour of similar behaviour for all L < N. While not constituting a fully conclusive proof of the 
converging overlaps, our results do demonstrate a striking similarity between the analytic structure 
of the exact ground state wave functions at L < N, and that of their CF counterparts. Results are 
^~j given for two different projection methods commonly used in the CF approach. 

PACS numbers: 03.75.Nt, 05.30.Jp, 71.10.Pm 



O 

(N 



I. INTRODUCTION 

The behaviour of ultracold atomic Bose gases under rotation has been a subject of intense study over the past years[T]. 
W) One of the most studied aspects, involving a large body of theoretical work, is the expected occurrence of strongly 

^-> correlated states of the fractional quantum Hall (FQH) type in the limit of ultrafast rotation[2 . On the experimental 

side, progress has been extremely impressive, from the creation of the first vortex in 1999 13] to arrays of hundreds of 
vortices[4]. Nevertheless there are practical obstacles to reaching the actual quantum Hall regime this way. Other, 
^h more promising, scenarios are now being pursued to try and create the effective 'magnetic' fields required for producing 
bosonic FQH states[5]. Another limit that has attracted considerable interest is that of slow rotation, in the sense 
that only one or a few vortices are present in the system. In particular, various groups have studied the yrast Zine[6 
of the system, i.e. the states with lowest energy for a given total angular momentum. Theoretical work includes 
analytical studies[B], exact diagonalization[7], the Gross-Pitaevskii approach [SJ [5], as well as trial many-body wave 
functions [ini E]. Moreover, for the lowest angular momentum regime, 2 < L < N, exact analytical ground state 
wave functions were found for bosons with short-range interaction sufficiently weak to assume the particles reside in 
the lowest Landau level [T2"Hro] . 

i— —i The focus of the present paper is this latter case. The composite fermion (CF) approach[T7], originally designed 

to give trial wave functions in the fractional quantum Hall regime, was first applied to some low-angular momentum 
states in Refs.|10[ lllj. It was notedpT] that overlaps between CF and exact wave functions for the 'single vortex' 
state at L = N appeared to increase with particle number, at least for small systems. This point was studied more 
systematically in Ref.[TH], and overlaps for the single vortex state computed numerically for up to 45 particles. The 

ly-j results strongly suggested convergence of the overlaps to unity in the thermodynamic limit, with deviation decreasing 
as ~ 1/N. This kind of behaviour is generically unusual for any kind of non-exact trial wave functions. In this case 

ly-^ it was even more surprising because a priori, the composite fermion construction could not be expected to work well 

f") in this low angular momentum regime - for reasons to be explained below. Unfortunately, the numerical calculations 
gave little insight into the reasons for this behaviour. Therefore, we here reexamine the issue analytically, including 

t-H also the states below the single vortex: We study the yrast states in the entire regime 2 < L < N, comparing the 

analytic form of the CF wave functions to their exact counterparts, and present evidence that the increase of overlaps 
with system size occurs for all these states. We address the issue for two different projection methods commonly used 
in the CF approach. In the case of 'full projection' (method I) we demonstrate that the analytic structure of the 
CF wave functions is strikingly similar to that of the exact ones; the latter have the form of fundamental symmetric 
polynomials in the coordinates (zi — Z), where Z is the center of mass, while the former are merely 'lacking' one or 
several coordinates in the center of mass term. For the other projection method (method II), it is possible to rewrite 
the CF wave functions, as well as the exact ones, in the form of an N x N determinant, divided by a Jastrow factor. 
We analytically compute overlaps between pairs of corresponding entries in these determinants and show that the 
overlaps between all pairs of entries converge to unity in the thermodynamic limit. This does not necessarily allow us 
to conclude that the overlaps between the full determinants (and thus the wave functions) go to unity as N — > oo, 
since the size of the matrices themselves grows with N; the latter remains a challenging open problem. 

Although exact analytic yrast wave functions are known for the particular states addressed in this paper, it is 
nevertheless of interest to examine why the CF trial wave functions work so surprisingly well in this regime. In 
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particular, this may help justify applying the CF construction to nearby states where no exact solutions are known - 
such as excitations above the present yrast states, or other few-vortex yrast states at L > N . 

The outline of the paper is as follows: Section [TT| presents some necessary background on rotating bosons in the 
lowest Landau level, the form of the exact yrast wave functions for 2 < L < N, as well as the composite fermion 



approach. Sections III and IV contain results for the single vortex state L — N for the two different projection 
methods. The methods developed here can be directly transferred to all the lower angular momentum states, L < N; 
this is the topic of section [V] Section |VI| gives some concluding remarks. An appendix is included, providing some 
details of the overlap calculations in section |IV| 



II. BACKGROUND 
A. Rotating bosons in the lowest Landau level 

Consider a system of N spinless bosons with mass to in a harmonic trap of strength to, rotating with angular frequency 
fi and interacting via a short-range (delta function) potential Hj . In a rotating frame the Hamiltonian can be written 
as 
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N v # i 

2m 2 



j=l 



ttL z + H T (1) 



where L z denotes the angular momentum around the rotation axis. Completing the square inside the brackets, this 
may be recast as 



N r 2 



+ (w - n)L z + Ht (2) 



with A = muj(—y,x), Hho{z) denoting the z-part of the harmonic oscillator Hamiltonian and || denoting the planar 
(x, y) part of the Hamiltonian. We note that the planar part of H takes the form of particles moving in an effective 
"magnetic" field B e ff =Vxi = 2mojz. 

Now, the interaction is assumed to be weak in the sense that it does not mix different harmonic oscillator levels. 
We will be interested, for a given total angular momentum, only in the lowest many-body states (the "yrast" states). 
In this limit, the model may be rewritten as a lowest Landau level (LLL) problem in the effective "magnetic" field 
B e ff = 2muj (and of course, n z = for the harmonic oscillator in the z-direction). The Hamiltonian then takes the 
form 

H = (oj-n)L + gY / 5 2 (ri-r i ) (3) 

i<j 

(?i = 1), where we now use L to denote the total angular momentum, L — J^ U = L z . 

In the following we shall be concerned with many-body wave functions describing this system at some given number 
of particles and total angular momentum. A convenient basis of single particle states spanning our Hilbert space (the 
lowest Landau level) is given by 

r) Ql = - l =z l e^- z/4 (4) 

where z = v 2moj(x + iy) are complex coordinates denoting the particle positions in the plane, and I is the angular 
momentum of the state. A general bosonic many-body wave function ip(zi, ...zjy) will then be a homogeneous, 
symmetric polynomial in the zt's, times the exponential factor exp(— ^2 i \zi\ 2 /4) (which will be suppressed throughout 
most of this paper for simplicity). The degree of the polynomial gives the total angular momentum of the state. 

In connection with the CF construction, we shall need more general single particle basis states for the full Landau 
problem, i.e. all Landau levels. In symmetric gauge, these are given by 

n =N e~ |z|2/4 z m £ m (— \ (5) 

where n is the Landau level index, to denotes the angular momentum, N nm is a normalization factor, and L™ are 
the associated Laguerre polynomials. 



B. Exact ground states 

Here and in the following we will use the notation ip[L,N] f° r the (ground state) wave function of an iV-particle state 
at total angular momentum L. Exact LLL ground state wave functions for the Hamiltonian discussed above are 
known for all angular momenta 2 < L < N. As was shown in a series of papers fT^HTB] . the yrast states in this 
angular momentum interval are simply given by fundamental symmetric polynomials Sl{z{) where Zi — Zi — Z, and 
Z = (J^i Zi)/N is the center-of-mass coordinate: 

^[L,N]= Yl ( z Pi~ Z )( Z P2- Z )---i z PL- Z )- (6) 

Pl<P2<-<PL 

For example, ip[2,N] = <5 [{z\ — Z)(z2 — Z)\, with S denoting symmetrization over all particle coordinates. As a special 
case the so-called single vortex state at L = N obeys the exact ground state wave function 

N 



C. Composite fermion approach 

The phenomenology of composite fermions\17\ was first introduced in the context of the fractional quantum Hall effect 
(FQHE), where it was used to construct very successful trial lowest Landau level many-body wave functions for a 
large number of FQH states. Roughly speaking, the basic idea of this construction is that an even number of vortices 
is bound to each electron. Each such vortex, in a mean field sense, cancels one flux quantum of the external magnetic 
field, so that the strongly interacting electrons are mapped to weakly interacting composite fermions, moving in a 
reduced effective magnetic field. A mathematical way of seeing that these composite objects are weakly interacting 
is that 'attaching p vortices' amounts to multiplying the many-body wave function by ni<i( z i — z j) P ~ a factor that 
'keeps the particles apart', or introduces an additional p-fold correlation hole around each particle, since it can be 
seen to go to zero as any two particle coordinates approach each other. A modified version of this model was later 
applied to study quantum Hall-type states in rapidly rotating Bose gases, mapping the lowest Landau level bosons to 
composite fermions by attaching an odd number of vortices [2l ITU1 fTTl ITUl |2"U] . Due to the above qualitative picture of 
composite fermions being weakly interacting objects, they are commonly approximated as non-interacting. Trial wave 
functions for LLL bosons at some angular momentum L are thus constructed as a Slater determinant $ of composite 
fermions with angular momentum L — pN(N — l)/2, times an odd power p of Jastrow factors, 

i®v] ({*»=? (*({*•*}) n(*i-* J -) p ) • (8) 

The operator V projects the wave function to the lowest Landau level. This is necessary since, in the intermediate 
step, the Slater determinant $ typically contains particles occupying states in higher Landau levels in the effective 
magnetic field, cfr. Eq.([5]). Essentially, the projection amounts to replacing z% by d/dzi in the polynomial part of 
the wave function. In the quantum Hall literature there are several ways of doing this in practice, and numerically 
the results are known to depend very little on the projection method chosen[T7]. We shall get back to the details 
of this in later sections. However, let us point to one important issue already now: Assume the Slater determinant 
<& is such that a majority of the composite fermions reside in higher effective Landau levels. This would translate 
to large powers in z (see Eq.([5])) and thus, upon projection, to a large number of derivatives acting on the factor 



Eli ( z * — Zj) p . Given the above qualitative arguments, if the number of derivatives is sufficiently large compared to 
the degree of the polynomial itself, one might expect that these derivatives would destroy the good correlations built in 
by the 'vortex attachment', so that the approximation of non- interacting composite fermions underlying ^ might be 
rather poor. In particular, naively applying the recipe (1SJ) to the lowest angular momentum states at L < TV involves 
0(N 2 ) derivatives acting on an order N 2 polynomial in Zi, killing off basically the entire vortex attachment factor, to 
leave a polynomial of degree N or less. So these states are far outside the quantum Hall regime, and one would not, 
a priori, expect the CF scheme to work very well. Nevertheless one finds the very surprising fact that the overlap of 
the composite fermion trial wave functions with the exact ones (k)]) increases with increasing particle number (being 
very large already for small systems), and approaches 1 as N — > oo - in other words, the composite fermion wave 
functions appear to be analytically exact in the thermodynamic limit. For L = N, this trend has previously been 
pointed out in the context of numerical calculations |11| [18]. In the present paper we study the analytic structure of 
the yrast wave functions in the entire interval 2 < L < N, and present analytic evidence of how this highly unusual 
behaviour comes about. 
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FIG. 1: Numerically computed overlap between exact and CF single vortex wave function for projection method I|18j. 

III. SINGLE VORTEX - PROJECTION I 

With the preliminaries in place, we are ready to study the single vortex state, L = N . For 'full projection' (referred 
to here as projection method I and explained below), this case was previously studied numerically [IS]; in Figjljwe see 
numerically computed overlaps between the composite fermion and exact single vortex wave functions for up to 40 
particles. This figure shows how the overlap converges towards unity, which is quite unusual behaviour for trial wave 
functions in general. As argued above, the fact that we are far outside the regime where the CF ansatz is naively 
expected to work, makes this result even more surprising. A full analytic proof for the convergence of these overlaps 
is still missing. But in this section we show how the CF wave function may be recast in a closed form, demonstrating 
that the analytic structure of the Jain states is extremely similar to that of the exact ground state; the main difference 
lies in that the center of mass coordinate in the exact wave function is replaced by 'incomplete' center of mass factors, 
with one or more particle coordinates 'missing'. The results of the present section can be directly used to find similar, 
closed expressions for the lower angular momentum Jain states, L < N. These results will be presented in section IV] 
According to Eq.((7]), the exact ground state wave function at L = N is given by Ili^i — %) where the center of mass 
coordinate Z — Y^Zi/N. The bosonic CF ansatz for the single vortex state, on the other hand, is given bv |lUlTS] 



1>[N,N\{{*i}) = V *[iV,JV] ({*.*}) II( Z * - z ^ ' 



(9) 



l<3 



where the unprojected CF Slater determinant is 
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Now, what we call projection method I, or full projection, simply amounts to replacing all z's in the Slater determinant 
by derivatives, z\ — > df, thus acting on the full Jastrow factor to its right. We now present an iterative method by 
which we can perform all the derivatives to find closed, analytic expressions for any number of particles. To this 
end, first note that for L = N = 3 we know from translation invariance arguments [TS] that the CF ansatz has to be 
identical to the exact wave function, 
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= [z x - Z) (z 2 - Z) (z 3 - Z) 



(11) 



with Z = (z 1 +z 2 +£3)/3, as can of course be checked explicitly. This will help us simplify the corresponding expressions 
for larger N. To see how this works, consider the case L = N = 4 and expand the 4x4 Slater determinant by the 



row with the largest number of derivatives, 

« = £(-ir +1 
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where the superscript (i) on the Slater determinant denotes that coordinate i is excluded. Next, note that the Jastrow 
factor in ( 12 ) may be split up according to 



4 (i) 3 

n (*«-*/»)= n (^-z P )-(-iy +i i[(z t ~z k ), 

a</3 a<P k^i 



(13) 



where the sign factor comes from moving coordinate i to the front. Furthermore, up to irrelevant overall multiplicative 
factors, 






(14) 



fc^i 



where Z is the center of mass of all four particles, while Z^ = Yljjti z j/3 is the three-particle center of mass with 
coordinate i omitted. Combining Eqs. (12), (13 1 and (14 1, we find 
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(15) 



The derivatives in the reduced Slater determinant might in principle act through the Jastrow factor, onto the factor 
(zi — Z). However, it can be shown that any such contribution vanishes due to antisymmetrization. Thus, we identify 
the first part of (15) as the three-particle single vortex wave function (11), giving 
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or, equivalently (again neglecting irrelevant overall multiplicative factors), 
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(17) 



Note that, except for the one missing coordinate in the center of mass, this is identical to the exact ground state wave 
function tf?\. If this was the case at all N, it would thus not be very surprising that the CF wave function approached 
the exact one in the thermodynamic limit. However, as we shall see now, going to higher particle numbers introduces 
factors where more and more coordinates are missing from the center of mass, as compared to the exact state. Thus, 
the convergence of the CF state to the exact one appears to be quite subtle. [23] 

The calculation outlined above for four particles, can be immediately iterated to larger systems. For example, for 
six particles, 
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It is straightforward to derive the general result 
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which can be expanded iteratively in analogy with Eq.(19|. Thus, the structure of the CF wave functions is remark- 



ably similar to the exact wave functions, the only difference being the incomplete center of mass factors. [25] While 
suggestive, this is not a proof that the CF wave functions approach the exact ones in the large N limit, since the 
number of coordinates 'missing' in the center of mass piece increases with N. Unfortunately, full analytic overlap cal- 
culations on these expressions seem quite intractable. Instead, we now revisit the problem within projection method 
II, where a partial analytic proof of the converging overlaps is possible. 



IV. SINGLE VORTEX - PROJECTION II 

In this section we reexamine the issue using a different LLL projection method, referred to as method II (see below 
for details). This is the most commonly used projection method in the CF quantum Hall literature, being numerically 
more manageable (considerably fewer derivatives) than full projection. For quantum Hall states, it is known to make 
basically no difference numerically which of these projections is used|17j. In our case, this is not entirely obvious; the 
increase and convergence of the overlaps is a subtle effect, and exactly how the derivatives act in the projection, may 
well make a difference. In order to compare the composite fermion ansatz to the exact wave function, we here rewrite 
each of them in the form of an N x N determinant divided by a single Jastrow factor. This has the advantage that 
we can, in fact, compare the two determinants entry by entry. We will prove analytically that for each pair of entries, 
the overlaps typically decrease with N for small numbers of particles, but then start increasing monotonically , and 
all converge to 1 in the thermodynamic limit. While looking very suggestive, this however does not necessarily prove 
that the same would be true for the entire determinant, i.e. for the overlaps between the full wave functions. This 
point will be discussed in more detail below. 

A. A useful identity 

We start this section by deriving a mathematical identity that will turn out very useful in the following. Consider the 
Slater determinant 
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which includes the Vandermonde determinant as the special case k = N. In general, being fully antisymmetric, A( ' 
is always divisible by a Vandermonde determinant, 



A^[z 1 ,...z N ] = f^_ k ({z i })-l[(z i -z j ), 



(22) 



l<3 



where /jy_/ £ is some fully symmetric polynomial of degree N — k. In fact, /iv-fc i s ; by definition, the Schur 
polynomial [2"T] s\({zi}) where the partition A = (I, I, ...I, 0, 0, ..., 0) contains (N — k) I's followed by k zeroes. From 
the Jacobi- Trudy identity [21]. which prescribes how to express Schur polynomials in terms of fundamental symmetric 
polynomials, it furthermore follows that f^_^({zi}) simply equals Siy-k({zi}), the fundamental symmetric polynomial 
of degree N — k in the N coordinates {zi}. Thus, 
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Some special cases of this identity have been discussed in the quantum Hall literature, e.g. 



Exact state 



Again, the exact ground state wave function at L — N is given by 5/v({,Si}) where z% = z% — Z. Using Eq.(23l with 
k = 0, and Zi replaced by Zi (and noting that (zi — Zj) = {zi — Zj)), we can write this wave function in the form 
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with 
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z\ z 2 
This form will turn out to be particularly convenient when comparing to the Jain wave function. 



C. Composite fermion ansatz 



The unprojected CF wave function for the single vortex was given in section III To put this into a form comparable 
to Eq.(24), we use a standard trick due to Jain|17j to absorb two Jastrow factors into the determinant and perform 
the projection z — > d entry by entry, 
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This is what we refer to as projection method II. 



D. Overlaps 
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To summarize the above manipulations, we have put both the exact and CF wave function into the form 
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where M is an N x N matrix. Recalling that our goal is to shed light on to what extent these two wave functions 
appear to approach each other with increasing N, we thus compare the two determinants, and do so entry by entry: We 
analytically compute pairwise overlaps (Mf^\M^ F ) between corresponding individual entries of the two determinants, 
and show that each one of these overlaps converges towards 1 in the limit N — > oo. In other words, the two matrices 
M ex and M approach each other entry by entry in the thermodynamic limit. Now, since in this limit the size 
of the matrices themselves goes to infinity, this does not strictly speaking prove that the overlap between the full 
determinants does approach unity as well. We have, so far, not been able to come up with a conclusive argument to 
settle this mathematically highly nontrivial question. 

In order to compute the overlap between LLL polynomials, one uses the orthonormality relation 
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for all coordinates. The actual overlap calculations involve a rather large amount of straightforward but tedious 
algebra and combinatorics. The main steps are outlined in appendix [XJ here we merely state and analyze the results. 



As a concrete example, let us start with the first row and compute the overlap between Aff^ and M^f as defined 
through (251 and (27), respectively. The result will only depend on the row index, i.e. be independent of the second 



index a, so we compute the overlap 
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where the denominator has to be included since the wave functions are not a priori normalized. The result is 
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Now, using the fact that lim„_ i . 00 (1 + l/n) n = e, and likewise, 5^^_ (l/fc!) = e, one finds that 
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So the overlap between any two corresponding entries in the first row approaches unity in the thermodynamic limit. 
Similarly, one finds for the second row, 
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which is again seen to converge towards unity in the limit N —> oo. In fact, as shown in the appendix, one can derive 
an analytic expression for the overlaps in any row a + 2, where a — 1, • • • , (A — 2), 
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We already examined some of the lowest rows. Let us also have a look at a couple of examples for the highest row 
indices. First, note that for the last (iVth) row, the overlap equals unity for all N. This can be shown by evaluating 
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Eq.(34| for a — N—2, or by noting through direct calculation that d 1 rife^i( z i — z k) = X)i=2( z i~ z «) = N(z\ — Z) a 



'or row number N — 1, i.e. a = N — 3, Eq.(34) gives, after some algebra, 
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which is again seen to approach unity as N — > oo. Similarly, 
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etc. Figures [2] and [3] show some examples of overlaps versus particle number N. They illustrate that O n +i > O 



N 


L 


Projection I 


Projection II 


4 


4 


0.980 


0.944 


5 


5 


0.986 


0.917 


6 


6 


0.990 


0.894 



TABLE I: Numerically computed overlaps between the full exact and CF wave functions at L = N for the two projection 
methods discussed in the text. Note that the overlaps increase with iV for full projection (method I), but, for these small 
particle numbers, appear to decrease with N for projection method II. 



for given N, i.e. overlaps increase with increasing row index, and the first row, Eq.(31), represents a lower bound 



of how fast the overlaps approach unity for growing A. As mentioned above, this does not provide a conclusive 
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FIG. 2: Some row overlaps, given in Eqs.(31 1- ( 34 1 , as functions of particle number. From below, rows 1 and 2 (too close to be 
resolved from each other), 3, 4, 5 and 6. 
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FIG. 3: Overlaps for first row (lower curve) and row number 600, for up to 4000 particles. 

answer concerning the overlaps between the full wave functions. In fact, from table IT] we note that at small particle 
numbers, projection method II produces a decreasing overlap between the full wave functions as function of N, which 
may look discouraging at first sight. However, as figure [2] illustrates, for system sizes this small, overlaps between 
the rows themselves actually tend to decrease rather rapidly as well, before passing a minimum and then increasing 
monotonically towards 1. It is thus fully conceivable that the overlaps between the full wave functions start increasing 
at larger N, too - but the present analysis does not provide a conclusive answer to this question. It would be 
interesting to study this point numerically for rather large systems, say, 20-30 particles or more. Incidentally, if it 
did turn out that the full overlaps do not converge towards unity at large N, this would be a very rare example of a 
system where the choice of projection method makes an important qualitative difference to the behaviour of the CF 
wave functions. 



V. LOWER ANGULAR MOMENTA 



The techniques developed for the single vortex in the previous sections, can be immediately transferred to the yrast 
states below the single vortex, L < N, and we shall again do this for the two projection methods separately. One 
apparent complication arising here is that for L < N , there will be more than one candidate CF Slater determinant for 
each state, the number of candidates increasing with decreasing L. So naively one would expect the ground state to be 
some linear superposition of these. However, remarkably, it turns out that after projection and antisymmetrization, 
the different CF candidates for a given L and N typically result in identical wave functions, leaving us with only one, 
representative, candidate state. 
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A. Projection I 



We start by discussing the case L = N — 1, before generalizing to lower angular momenta. To construct this state 
using the CF scheme, one finds two candidate Slater determinants for the ground state|27j: 
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Manipulating these, one can show that these two candidates, when acting on a single Jastrow factor, produce identical 
wave functions. [28. Thus, it is sufficient to consider one of them, and we shall choose the first, Eq.(37). Proceeding 
similarly to the steps of Eqs. (12 1 - (16), we expand the Slater determinant (37) by the row with the highest number 



of derivatives to cast the CF wave function in the form 
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(39) 



Noting that the A^ — 1 derivatives simply turn the last part into a constant, and that the first part has the form of a 
sum over (N — l)-particle single vortex states (discussed in section III) with particle i omitted, we thus have 



iV 



ih CF 

V[N- 



1,N] 



/ , r[jv-i,jv-i]- 



(40) 



So the CF yrast state at L = N — 1 is simply a symmetrized sum over L = N states with one particle less; these, in 
turn, are known explicitly from section \UJ\ 

Similar considerations apply for lower L, where there will be even more CF candidates for each state. Although we 
have not bothered to perform a complete, systematic study for all possible cases, we have checked for a large class of 
these candidates that they again produce identical wave functions (or zero) . Since the analytic structures of the wave 
functions at lower L are very analogous to those encountered above, we believe this to be true in general 29J. We thus, 
again, concentrate on one representative candidate, most con veniently the one whose Slater determinant has rows 
[z, 1, d, ...,d N ~ k ~ 2 ,d N ~ k , ..., c^^ 1 ] , the generalization of (37). The above calculation then immediately generalizes 
to give the result 
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[N-k,N-l]- 



(41) 



Thus, using this expression iteratively, the CF yrast state at any L — N — k can be expressed as a symmetrized 
(multiple) sum over single vortex states with N — k particles, which are known from section III For example, 
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and analogously for smaller L. Note again the striking similarity to the structure of the corresponding exact yrast 
wave function ([6]) . Table fil] shows some numerically computed overlaps between the exact and yrast wave functions 
for the L = N, N — 1, and N — 2 states. We note that in all cases, the overlaps indeed increase with the number of 
particles, just as for the single vortex state[3D]. At the same time, overlaps tend to decrease as one moves from the 
single vortex towards lower angular momentum states at fixed L or TV. This is perhaps not very surprising, given that 
the state at L = N — k can be expressed through single vortex states at a smaller number (N — k) of particles. The 
results of this section thus strongly suggest that the same kind of behaviour will be found for all yrast states in the 
interval 4 < L < N (not including L — 2,3 here since there the CF wave function is known to be exact for all iV[ll]). 



L 


L=N 


L=N-1 


L=N-2 


4 


0.980 


0.953 


0.926 


5 


0.986 


0.967 


0.945 


6 


0.990 


0.973 


0.954 



TABLE II: Numerically computed overlaps (projection I) between exact and CF wave functions at L — N, N — 1, and N — 2 for 
different values of L (so the values of N range from 4-6, 5-7, and 6-8, respectively in the three columns). Note that the overlaps 
increase with N in all three cases. For fixed L, they tend to decrease as one moves away from the single vortex towards lower 
angular momenta. 



B. Projection II 



The analysis of the L < N yrast states with projection II will be very analogous to that of section |IV| for the single 
vortex. The aim of this section is to demonstrate how we can, again, express both the exact and CF states as a 
determinant divided by Jastrow factor, and compare the determinants entry by entry. Considering first the exact 
ground state at L = N — k, Eq.([6|, we can again use Eq.(23), with general k, and Zi replaced by Zi, to write 
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For the CF yrast state at L = N — k, we again choose to study the representative Slater determinant discussed in 
connection with Eq.(41 1. This immediately leads to the following analog of Eq.(p7|, 
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Note that the last row of M&f_ k N , reduces to a row of identical constants, thus being identical (up to an overall 
constant) to the last row of Mffi_ k N ,. All other rows have entries we already encountered in our discussion of the 
single vortex in section [rv] From the results of that section, we can thus immediately conclude that between each pair 
of entries of these two determinants, for all 4 < L < N , the overlaps converge to unity in the thermodynamic limit. 
In other words, all yrast states below the single vortex behave in a way very analogous to that found for L = N . 

VI. CONCLUSIONS 

The aim of this paper was to shed light on the surprising fact that trial CF wave functions for yrast states at and below 
the single vortex increase with system size, and appear to converge to the exact ground states in the thermodynamic 
limit, as suggested by numerics using projection method I. The ultimate goal would of course be a full derivation in 
terms of analytic expressions for the overlaps between CF and exact wave functions for arbitrary N . One might also 
wish to study analytic expressions for the eigenenergies of the CF states as functions of L and N and compare to 
their exact counterparts[16 . Both tasks turn out to be very challenging. Here, we reported progress in this direction 
by comparing analytically the form of the CF and exact yrast wave functions for the two most common projection 
methods, and pointing out their striking similarities. In particular, our analytical overlap calculations (in projection 
method II) for individual matrix entries should be viewed as suggestive evidence of the converging overlaps of the 
full wave functions, but do not provide a conclusive answer - the two main 'problems' being that the sizes of the 
matrices at hand themselves go to infinity in the thermodynamic limit, and that the increase of the individual row 
overlaps is non-monotonic. It would thus be interesting to study the full overlaps for projection method II numerically 
for large systems. If they do not converge to unity, this would be one of very few known cases where the choice of 
projection method in the CF construction makes an important qualitative difference. One may also hope that some 
of the identities and manipulations pointed out in this paper, may turn out useful in other contexts, where exact 
analytic wave functions are not available. Finally the calculations in this paper were done entirely in disk geometry; 
it would be of interest to reexamine the issues discussed here in, e.g., spherical geometry. 

Acknowledgement: We thank Stellan Ostlund, Steve Simon, S0ren Eilers, Hans Hansson, and Bcrnhard Mchlig for 
enlightening discussions. This work was supported by the Norwegian Research Council. 

Appendix A: Analytic overlap calculations 



In this appendix, we derive the expression for the overlap given by Eq. ( 34 ) for a general number of particles N and a 



given row j3. To prove this result we extensively use the binomial and multinomial expressions : 

(a - b) a = ^(-1) M ( a J a l b°" 1 (Al) 
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\ i / |k|=a \ / i 
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IW 



where k = (k\, ■ ■ ■ , k n ) and |k| = ^" fej = a. Since we are interested in the overlap between the matrix entries MS* 
and M£ F , defined through Eqs. (25 1 and (271 respectively, we first calculate the norm of each polynomial. The norm 
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of MS* is given by (consistently suppressing factors of I-k from here on) 
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In the first step, we applied the above binomial expansion and the orthogonality relation ( 29 1 for the variables z%. Then, 



Eq. ( A7 1 can be obtained after expanding the remaining sum with the multinomial expression and the orthogonality 



relation. Further simplifications give rise to Eq. (A9). 

To calculate the norm of M£ F , we first need to evaluate the following expression: 
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where S a ({zi}) are the elementary symmetric polynomials of degree a. The norm of Eq. (AlOl can then easily be 
deduced by using the orthogonality relation p9|. It is given by 
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Finally the overlap is given by 
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To obtain Eq. (A15l, we first applied the orthogonality relation to the variable Z\, then expanded the remaining 



multinomial expression, only keeping terms involving k a = or 1 (as other terms do not give any contributions to 
the overlap). Combining Eq. (A17), (A13) and (A9|, we thus find that the correctly normalized overlap for any given 
row (3 + 2 is given by 
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By reorganizing the factorials and products in the denominator of Eq.( A18), the sum over a can be recast to give the 
final expression (34). 



The overlaps 0\ (31) and O2 ([33]) for the first and the second rows (where the CF determinant does not contain 
derivatives) are calculated separately, but the calculations follow the same general strategy as those demonstrated 
here. 
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